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Oa ' The interstellar medium is structured as a hierachy of gas clouds, that looks self- 

similar over 6 orders of magnitude in scales and 9 in masses. This is one of the more 
extended fractal in the Universe. At even larger scales, the ensemble of galaxies 
looks also self-similar over a certain ranges of scales, but more limited, may be over 
3-4 orders of magnitude in scales. These two fractals appear to be characterized 
by similar Hausdorff dimensions, between 1.6 and 2. The various interpretations of 
O^ ' these structures are discussed, in particular formation theories based on turbulence 

^Nj ' and self-gravity. In the latter, the fractal ensembles are considered in a critical 

state, as in second order phase transitions, when large density fluctuations are 
observed, that also obey scaling laws, and look self-similar over an extended range. 
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1 Introduction 

Fractals are ensembles that can be defined by their self-similarity. The name 
has been introduced by Mandelbrot (1975) to define geometrical or mathemat- 
ical sets, that have a non-integer, i.e. fractional dimension. The dimension 
determines whether a system is homogeneous, and what fraction of space is 
filled. For a homogeneous density, the mass of the medium is increasing as 



j^ , the 3rd power of its radius (in 3D) , while when the medium is fractal, he may 

occupy a tiny fraction of space, and the mass contained within a scale r is 

- ^-, , Moc r^ , with D, the Hausdorff dimension, lower than 3. 

^^ ' The two fractals described here, the interstellar medium (ISM) and galax- 

^ , ies, have a dimension D w 1.7. Of course these physical ensembles are only 

approximations of mathematical fractals. They are self-similar only between 

two limiting scales, where boundary effects occur, while a pure mathematical 

fractal is infinite; and they arc quite randomly distributed, their self-similarity 

being only statistical. 
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Figure 1. Left: IRAS 100 fj,m map of the Taurus molecular cloud complex, traced by the 
dust emission. The square is ~ 4000 pc'^ . Right Zoom of the central region (the square is 
now ~ 400 pc^). 



2 The Interstellar Medium 

The gaseous interstellar medium has a very irregular and fragmented struc- 
ture. It consists of clouds of hydrogen, either atomic or molecular, according 
to its density or column density. The atomic gas is less dense and more diffuse 
in general, while the molecular gas gathers the most clumpy, cold and dense 
phase (see the survey of the Milky Way in Molecular Clouds by Dame et al 
1996). 

At any scale, the self-similar appearance of the clouds ensemble makes it 
difficult to determine their absolute scale on photographs, without any other 
information (about velocity, distance, etc.). In Fig |l|, we show for example 
maps of the nearby Taurus cloud, observed in the far-infrared at lOO/xm with 
the IRAS satellite. The emission comes from dust heated by nearby stars, 
or by the interstellar radiation field. The right-hand side map is an enlarged 
view of the central left map, but still quite similar. Constraints due to spatial 
resolution of telescopes, cannot in general allow to observe more than a few 
orders of magnitude in scale range for the same cloud, but the observation of 
several clouds at various distances in the Milky Way (from 0.05 to 15 kpc) or 
in external galaxies, have established the scaling laws over a wider range. 
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2.1 Low and High-mass Cut-off for the ISM Fractal 

The largest self-gravitating entities in the Galaxy are the so-called Giant 
Molecular Glouds (or GMG) of about lOOpc diameter, and 10® Mq in mass. 
Larger clouds cannot exist since they would be teared off by the galactic shear, 
i.e. the tidal forces due to the galactic potential itself. This is the high cut-off 
scale in the fractal structure. What is the smallest size? 

It is difficult to observe directly in emission the smallest structure, due to 
lack of spatial resolution and sensitivity. But structures of about 10-20 AU in 
size (i.e. ~ 10~^ pc) have been observed for a long time through scattering 
of the quasar light (Fiedler et al. 1987, 1994, Fey et al. 1996): clumps in the 
electronic density diffract the light rays from remote quasars, and produce 
an "extreme scattering event" (ESE) lasting for a few months, in their rapid 
motion (100-200km/s) just in front of the quasar. QSOs monitoring during 
several years has determined that the number of scattering structures is 10^ 
times as numerous as stars in the Galaxy. The problem of stability and life- 
time of these structures, with much higher pressure than surroundings, can be 
solved if they are self-gravitating (Walker & Wardle 1998); they are then of 
10"'^ M0 in mass, and have a gas density around 10^" cm~^. They correspond 
to the smallest fragments predicted theoretically (Pfenniger & Combes 1994). 
These structures are now observed in a large number directly, through VLBI 
in the vicinity of the Sun, through HI absorption in front of quasars (e.g. 
Diamond e tal. 1989, Davis et al 1996, Faison et al 1998). If this 10-20 AU 
size is adopted for the low cut-off scale of the fractal structure, the latter 
ranges over 6 orders of magnitude in size, and about 10 in masses. 



2.2 Scaling Laws 

To better quantify the self-similar structure, several works have revealed that 
the interstellar clouds (either molecular or atomic) obey power-law relations 
between size, linewidth and mass (cf Larson 1981). These scaling relations are 
observed whatever the tracer. The original one is the size-linewidth relation 
(directly derived quantities) , while the mass is only a secondary quantity, very 
uncertain to obtain, since there is no good universal tracer. The H2 molecule 
does not radiate in the cold conditions of the bulk of the ISM (10-15 K), 
since it is symmetric, with no dipolc moment. The first tracer is the most 
abundant molecule CO (10^* in number with respect to H2), but it is most of 
the time optically thick, or photo-dissociated. The relation between the sizes 
R and the line-widths or velocity dispersion a, can be expressed through the 
power-law: 
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with q between 0.3 and 0.5 (e.g. Larson 1981, Scalo 1985, Solomon et al 1987, 
cf Fig. g). There are some hints that molecular clouds are virialised (at least 
at large scale, since the masses are even more uncertain at low scales and high 
densities). If the virial is assumed at all scales, then: 

cr^ ex M/R 
and the size-mass relation follows: 

with D the Hausdorff fractal dimension between 1.6 and 2. It can be deduced 
also that the mean density over a given scale R decreases as 1/R", where a 
is between 1 and 1.4. 

Other interpretations are possible (see e.g. Falgarone 1998). Recently, 
Heithausen et al (1998) have extended the size-linewidth and size-mass re- 
lations down to Jupiter masses; their mass-size relation is Mcx r^'^^, much 
steeper than previous studies, but the estimation of masses at small scales is 
quite uncertain (in particular the conversion factor between CO and H2 mass 
could be higher, and the mass at small-scale underestimated). It appears also 
that the self-similarity of structures is broken in regions of star formation. 
The break is observed as a change of slope in mass spectrum power laws at 
about 0.05pc in scale in the Taurus cloud, corresponding to the local Jeans 
length (Larson 1995). In other regions, cither this break does not occur, or is 
occuring at larger scales (Blitz & Williams 1997, Goodman et al. 1998). 

Dense gas is molecular, and cold H2 molecules are difficult to trace 
(Combes & Pfenniger 1997). The tracer molecules such as CO are either 
optically thick, or not thermally excited (in low-density regions), photo- 
dissociated near ionizing stars, or depleted onto grains. Isotopic species, like 
^^CO or C^^O, are poor tracers also, because of selective photo-dissociation. 
The large range of scales is also a source of bias, and of overestimates of the 
fractal dimension: the small scales are not resolved, and observed maps are 
smoothed out. This process, which makes the fractal look as a more diffuse 
medium of larger fractal dimension, can also lead to underestimates of the 
mass by factors more than 10 (e.g. simulations in Pfenniger & Combes 1994). 

The observations of molecular clouds reveal that the structure is highly 
hierarchical, smaller clumps being embedded within the larger ones. This 
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Figure 2. Top: Size-linewidth relation taken from various sources: Dam86: Dame et al. 
(1986); Sol87: Solomon et al. (1987); Hcit98: Hcitliausen et al. (1998); MBM85: Magnani 
ct al. (1985); Will94: Williams et al. (1994); Fal92; Falgarone et al. (1992); W95: Wang 
ct al. (1995); Ward94: Ward-Thompson et al. (1994); Lem95: Lemme et al. (1995). An 
indicative line of slope 0.5 is drawn. Bottom: Mass-size relation deduced from the previous 
one, assuming that the structures are virialised. The line drawn has a slope of 2. 



structure must be reminiscent of the formation mechanism, through recursive 
Jeans instabihty for instance. Since we have no real 3D picture, it is however 
difficult to ascertain a complete hierarchy, or to determine the importance of 
isolated clumps and/or a diffuse intercloud medium. An indicator of the 3rd 
dimension is the observed radial velocity, which is turbulent without system- 
atic pattern. It has been possible, however, to build a tree structure where 
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each clump has a parent for instance for clouds in Taurus (Houlahan & Scalo 
1992). 



2.3 2D-Projection of the Fractal 

The projection of a fractal of dimension D may not be a fractal, but if it is one 
with dimension Dp it is impossible a priori to deduce its fractal dimension, 
except that 

Dp = D if D<2 

Dp = 2 if D > 2 (Falconer 1990). 

In 2D, the fractal dimension can be measured by computing the surface 
versus the perimeter of a given structure. This method has been used in 
observed 2D maps, like the IRAS continuum flux, or the extinctions maps of 
the sky. In all cases, this method converge towards the same fractal dimension 
D2. For a curve of fractal dimension D2 in a plane, the perimeter P and area 
A are related by 

P oc A^^/2 

Falgarone et al (1991) find a dimension D2 — 1.36 for CO contours both at 
very large (degrees) and very small scales (arcmin), and the same is found 
for IRAS 100/x contours in many circumstances (e.g. Bazell & Desert 1988). 
Comparable dimensions {D2 between 1.3 and 1.5) are found with any tracer, 
for instance HI clouds (Vogelaar & Wakker 1994). 



Turbulence 



3. 1 Theory 

The interstellar medium is highly turbulent. To quantify, we can try to es- 
timate the Reynolds number i?e, that separates the laminar (low Rg) from 
turbulent regimes: 

Re — vljv 

where v is the velocity, / a typical dimension, and v the kinematic viscos- 
ity. Turbulent regimes are characterized by R^ >> 10"^, implying that the 
advection term v . V v dominates the viscous term in the fluid equation. 
The turbulent state is characterised by unpredictable fluctuations in density 
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and pressure, and a cascade of whirls. In the ISM, the viscosity can be es- 
timated from the product of the macroturbulent velocity (or dispersion) and 
the mean-free-path of cloud-cloud collisions (since the molecular viscosity is 
negligible). But then the Reynolds number is huge (w 10^), and the presence 
of turbulence is not a surprise. 

This fact has encouraged many interpretations of the ISM structure in 
terms of what we know from incompressible turbulence. In particular, the 
Larson relations have been found as a sign of the Kolmogorov cascade (Kol- 
mogorov 1941). In this picture, energy is dissipated into heat only at the lower 
scales, while it is injected only at large scale, and transferred all along the hi- 
erarchy of scales. Writing that the energy transfer rate v^ /{r/v) is constant 
gives the relation 

V ex r^/^ 

which is close to the observed scaling law, at least for the smallest cores (Myers 
1983). The source of energy at large scale could be the differential galactic 
rotation and shear (Fleck 1981). This idealized view has been debated (e.g. 
Scalo 1987): it is not obvious that the energy cascades down without any 
dissipation in route (or injection), given the large-scale shocks, flows, winds, 
etc... observed in the ISM. The medium is highly supersonic (with Mach 
numbers larger then 10 in general), and therefore very dissipative. Energy 
could be provided at intermediate scale by stellar formation (bipolar flows, 
stellar winds, ionization fronts, supcrnovae..). Also, the interstellar medium is 
highly compressible, and its behaviour could be quite different from ordinary 
liquids in laboratory. However a modified notion of cascade could still be 
applied, leading to a Burgers spectrum, with 

w oc r ' 

(see e.g. Vazquez- Semadeni et al. 1999). 

The role of the magnetic field is still unclear. The intensity of the field 
has been measured through Zeeman line effects to be around a few /iG in 
dense clouds (Troland et al. 1996, Crutcher et al 1999), and the observations 
are compatible with the hypothesis of equipartition between magnetic and 
kinetic energy. The magnetic field therefore plays an important role in energy 
exchange between the various degrees of freedom, but cannot prevent grav- 
itational instabilities and cloud collapse, parallel to the field lines; although 
Gammie & Ostriker (1996) had found that magnetic waves could indeed pro- 
vide some support in low dimensionality, in more realistic simulations there is 
almost no difference between B = Q and strongly magnetized models, as far 
as the dissipation and energy decay rates are concerned (Stone et al. 1998). 
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Besides, many features of ordinary turbulence are present in the ISM. 
For instance, Falgarone et al (1991) have pointed out that the existence of 
non-gaussian wings in molecular line profiles might be the signature of the 
intermittency of the velocity field in turbulent flows. More precisely, the ^"^CO 
average velocity profiles have often nearly exponential tails, as shown by the 
velocity derivatives in experiments of incompressible turbulence (Miesch & 
Scalo 1995). Comparisons with simulations of compressible gas give similar 
results (Falgarone et al 1994). Also the curves obtained through 2D slicing of 
turbulent flows have the same fractal properties as the 2D projected images of 
the ISM; their fractal dimension D2 obtained from the perimeter-area relation 
is also 1.36 (Sreenivasan & Meneveau 1986). 

More essential, the ISM is governed by strong fluctuations in density and 
velocity. It appears chaotic, since it obeys highly non-linear hydrodynamic 
equations, and there is coupling of phenomena at all scales. This is also related 
to the sensibility to initial conditions that defines a chaotic system. The 
chaos is not synonymous of random disorder, there is a remarquale ordering, 
which is reflected in the scaling laws. The self-similarity over several orders of 
magnitude in scale and mass means also that the correlation functions behave 
as power-laws, and that there is no finite correlation length. This characterizes 
critical media, experiencing a second order phase transition for example. 



3.2 Turbulence Simulations 

A large number of hydrodynamical simulations have been run, in order to 
reproduce the hierachical density structure of the interstellar medium. How- 
ever, these are not yet conclusive, since the dynamical range available is still 
restricted, due to huge computational requirements. To gain in spatial reso- 
lution, 2D or even less (because of symmetries) computations are performed, 
but often the results cannot be generalised to 3D. 

It has been argued that self-similar statistics alone can generate the ob- 
served structure of the ISM, in pressureless turbulent flows without self-gravity 
(Vazquez-Semadcni 1994); however, only three levels of hierarchical nesting 
can be traced. When heating and cooling processes are included, and since 
the corresponding time-scales are faster than the dynamical tine-scales, the 
gas can be derived by a polytropic equation of state, as P ex p'', with 7 be- 
ing the effective polytropic index. The isothermal value 7 = 1 is one of the 
possibilities, between and 2 found in simulations. 

From the size-linewidth relation a ex i?^'^, and the second observed scal- 
ing law p ex R^^, it can be deduced that 
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and therefore, if the turbulent pressure P is defined as usual by 

dP/dp = cr^ 

it follows that 

P ex logp 

which is the logatropic equation of state, or " logatrope" . This behaviour has 
been tested in simulations (e.g. Vazquez-Semadeni et al 1998), but the logat- 
rope has not been found adequate to represent dynamical processes occuring 
in the ISM (either hydro, or magnetic 2D simulations). The equation of state 
of the gas would be more similar to a polytrope of index 7 « 2. But the 
results could depend whether the clouds are in approximate equilibrium or 
not (cf McLaughhn & Pudritz 1996). 

Vazquez-Semadeni et al (1997) have searched for Larson relations in the 
results of 2D self-gravitating hydro (and MHD) simulations of turbulent ISM: 
they do not find clear relations, but instead a large range of sizes at a given 
density, and a large range of column densities; the Larson relation is more 
the upper envelope of the region occupied by simulated points in the p - R 
diagram. They suggest that the observational results could be artefacts or 
selection efi^ects (existence of a threshold in column density for UV-shielding 
for example). 

MHD simulations have also explored the correlations between magnetic 
field strength and density, and field directions and filamentary morphology 
of clouds. Again, there is no clear correlation, only a larger scatter, with an 
upper envelope such as B varying as p^'^ (Padoan & Nordlund 1998, Padoan 
et al. 1998). 

Note also that chemical reactions network, combined with turbulence, can 
be the source of considerable chaos (Rousseau et al. 1998). 

4 Self-Gravity 

4-.1 Theory 

Although the ISM is a self-organizing, multi-scale medium, comparable to 
what is found in laboratory turbulence, there are very special particularities 
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that are not seen but in astrophysics. Self-gravity is a dominant, while it has 
not to be considered in atmospheric clouds for instance. It has been recognized 
by Larson (1981) and by many others that at each scale the kinetic energy 
associated with the linewidths balances the gravitational energy: clouds are 
virialized approximately, given their very irregular geometry. 

Self-gravitating gas in an isothermal regime (which is a good approxi- 
mation for the ISM), is known to be subject to instabilities. If we consider 
a sphere of gas confined in a box, and in contact with a thermostat, it will 
tend to follow an isothermal sphere, if the gas is hot enough. Below a certain 
critical temperature, there is no equilibrium any more, and the gas heats up 
when being cooled down, it is the gravothermal catastrophy, caused by nega- 
tive specific heat (Lynden-Bell & Wood 1968). Small sub-condensations could 
form, and the physics will be more complex, since in the asymptotic case of an 
isolated clump, it will evaporate in a large number of dynamical times. The 
environment is quite important, and the fact that the clumps can exchange 
mass and energy with surroundings as well (e.g. Padmanabhan 1990). 

Gravitational instability and cloud collapse is accompanied by fragmen- 
tation in a system with very efficient cooling, and this process can provide the 
turbulent motions observed. The theory was first proposed by Hoyle (1953) 
who showed that the isothermal collapse of a cloud led to recursive fragmen- 
tation, since the Jeans length decreases faster than the cloud radius. Rees 
(1976) has determined the size of the smallest fragments, when they become 
opaque to their own radiation. They correspond roughly to the smallest scales 
observed in the ISM (sizes of 10 AU, and masses of 10"'^ M©, see the physical 
parameters of the " clumpuscules" in Pfenniger & Combes 1994). 

The hierarchical fractal structure yielded by recursive fragmentation can 
be simulated schematically, in order to estimate the resulting filling factor, 
and the biases introduced by observing with a limited spatial resolution (e.g. 
Pfenniger & Combes 1994). Clouds are distributed according to the radial 
distribution of an isothermal sphere in r~^, and fragments also have positions 
selected randomly according to this radial distribution, and so on, hierarchi- 
cally. The number TV of fragments at each level, and the fractal dimension D 
chosen suffice to determine the size ratio between to imbricated spheres, i.e. 
— N'^/D , A sample of these fractals is shown in Fig. 0, for two dimensions 
2.2 and 1.6. 

Many other physical processes play a role in the turbulent ISM, as for 
instance rotation and magnetic fields. But they cannot be identified as the 
motor and the origin of the structure. Galactic rotation certainly injects en- 
ergy at the largest scales, but angular momentum cannot cascade down the 
hierachy of clouds; indeed if the rotational velocity is too high, the structure 
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Figure 3. Left: Simulations of a fractal of dimension D = 2.2 built by recursive fragmenta- 
tion; Right: Same for a fractal of dimension D = 1.6, the number of points are similar in 
the two cases. 



is unstable to clump formation (cf Toomre criterium, 1964), and the non- 
axisymmetry evacuates angular momentum outside the structure. Magnetic 
fields are certainly enhanced by the turbulent motions, and could reach a cer- 
tain degree of global equipartition with gravitational and kinetic energies in 
tbe virialised clouds. But they cannot be alone at the origin of the hierar- 
chical structure, gravity has to trigger the collapse first. Besides, there is no 
observational evidence of the gas collapse along the field lines, polarisation 
measurements give contradictory results for the field orientation with respect 
to the gas filaments. Therefore, although rotation, turbulence, magnetic fields 
play an important role in the ISM, they are more likely to be consequences of 
the formation of the structure. 



4-2 Self- gravitating Simulations 

Simulations of self-gravitating gas are very demanding, since large gradients 
rapidly set up, and spatial resolution must be adapted. A general rule is that 
the resolution is well below the Jeans length at any point, but the Jeans length 
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shrinks along the collapse. Artificial fragmentation can sometimes happen due 
to artifacts (see e.g. Truelove et al. 1997). To compensate for the limited 
spatial range, periodic boundary conditions are used, to simulate the ISM. 
Klessen (1997) and Klessen et al. (1998) have considered the fragmentation 
of molecular clouds, in 3D, starting with an homogenous cloud with small 
primordial fluctuations. These initial conditions are very similar to what is 
used in a cosmological background. The fluctuations are a gaussian density 
field with power spectrum P(k) ex k~^. After one free-fall time, the gas has 
evolved into a system of filaments and knots, some of them contain collapsing 
cores. To avoid the problem of spatial resolution, the condensed cores are 
then replaced by sink particles, simulating therefore a low cut-off scale. 

Klessen et al. (1998) compute the mass spectrum of clumps, which looks 
very similar to the observed one dN/dm oc m~^-^ in the ISM, at least over 
1.5 order of magnitude. The results however, depends still a lot on the initial 
conditions, and the power spectrum used (Semelin & Combes 1999). One of 
the problem of these simulations, in comparison with the ISM, is the absence of 
energy re-injection, and systematic motions. In a Galaxy, differential rotation 
and shear should provide both. 

Taking into Account the Shear 

One way to re-inject energy at large scale is the differential rotation of 
the Galaxy, as we have already emphasized. Computations of self-gravitating 
shearing sheets have already given interesting large-scale structures, that 
might be similar ro fractals. Toomre (1990) carried out self-gravitating sim- 
ulations on gas particles, that were regularly cooled, in a shearing environ- 
ment, and a quasi-periodical boundary conditions. In fact, the periodicity 
has to take into account the shear motions, and proper sliding of the sheet at 
larger and smaller radii is necessary in order to ensure continuity. The main 
result is a steady wave pattern, that sets quickly in through gravitational 
instabilities, and differential rotation, and is continuously renewed. Toomre 
(1990) was surprised himself of the efficiency of the mechanism, and of its 
quasi-stationarity: the energy dissipated in the gas cooling was compensated 
exactly by the rotational shear. This could be the main energy source was 
the gravitational structures in the ISM of galaxies. Toomre & Kalnajs (1994) 
refined the computations, and obtained the "typical" morphology of the pat- 
tern. These simulations demonstrate the power of gravity, allied with shear, 
to produce spiral structure, by the swing mechanism. They also show how 
much structures at all scales, or spiral chaos, can be sustained and maintained 
by the self-gravitation of a cooled granular distribution. 

The work was taken up recently by Huber & Pfcnniger (1999), that gener- 
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alised the computations in 3D, taking into account the galaxy plane thickness 
of the gaseous component. Their cooling is simulated by a viscous force, pro- 
portional to the particle velocity. They also find the characteristic filaments 
due to the shear combined with self-gravity, and measure the corresponding 
fractal dimension of the structures. Since the filaments are much longer in the 
plane than thick (perpendicular to it), the fractal dimension is in fact domi- 
nated by the z-dimension plane morphology, at least at large scales. So the 
range of scale where the fractal dimension is lower than 2 is very small. The 
boarder effects are then too large to determine D without ambiguity. These 
calculations are however encouraging that a fractal structure could develop at 
sub-kpc scales, when shearing is dominant, near the high cut-off of the ISM 
fractal. 

In a recent work (Semelin & Combes, 1999), we have also tried to de- 
termine the fractal dimension of a nearly isothermal self-gravitating gas. A 
tree-code and cloud-cloud collisions algorithm was used in 2D and 3D to fol- 
low the collapse of a periodically-replicated piece of ISM. Without shear, the 
collapse of initial perturbations are followed , and a fractal structure is found 
only in a transient way. With energy re-injection, and in particular with shear 
(and Coriolis forces), the structures are formed over several orders of magni- 
tude. The characteristic spiral filaments are observed, and embedded clumps 
form, and can disappear and reform stochastically. Simulations of a shearing 
sheet are shown in Fig |[ The fractal dimension derived is around 1.8. 

Wada & Norman (1999) have also simulated a shearing galaxy disk, but 
with a multi-phase medium: they obtain the formation of dense clumps and 
filaments, surrounded by a hot diffuse medium. The quasi-stationary filamen- 
tary structure of the cold gas is another manifestation of the combination of 
shear and self-gravity in a cooled gas. 



5 Galaxy Distributions 

It has been known for a long time that galaxies are not distributed homoge- 
neously in the sky, but they follow a hierarchical structure: galaxies gather 
in groups, that are embedded in clusters, then in superclusters, and so on 
(Charher 1908, 1922, Shapley 1934, AbeU 1958). Moreover, galaxies and 
clusters appear to obey scaling properties, such as the power-law of the two 
point-correlation function: 



^(r) ex r ^ 
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Figure 4. From Left to Right, and Top to Bottom: Simulations in 2D of a shearing sheet 
of gas particles, subject to their own self-gravity, and to collisions (dissipation). The central 
square in each frame is the area simulated, and the lanes surrounding it have been plotted 
to show the particular periodicity, that takes into account the shear (relative X-motions). 
The gradient of the rotational angular velocity Q is along the Y-axis. The time is 1, 4, 11, 
and 16 free-fall times respectively for the four successive frames. 



with the slope 7, the same for galaxies and clusters, of w 1.7 (e.g. Peebles, 
1980, 1993). According to the self-similar morphology, and the scaling laws, 
the galaxy ensemble can also be characterized by a fractal. But what are the 
limiting scales? 
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5.1 Low and High-mass cut-off for the galaxies fractal 

The smallest structure is of course a galaxy, and the starting point of the 
fractal is therefore obvious: 10 kpc in scale, and about 10"'^'' M©, representative 
of a dwarf galaxy. The high-mass cut-off is much less obvious, and there has 
been a continuous debate in the recent years to know what is the scale of 
transition to homogeneity, and even whether this scale exists. 

Isotropy and homogeneity are expected at very large scales from the Cos- 
mological Principle (e.g. Peebles 1993). The main observational evidence 
in favor of the Cosmological Principle is the remarkable isotropy of the cos- 
mic background radiation (e.g. Smoot et al 1992), that provides information 
about the Universe at the matter/radiation decoupling. At very large scales, 
the Universe must then be homogeneous. There must exist a transition be- 
tween the small-scale fractality to large-scale homogeneity. This transition 
is certainly smooth, and might correspond to the transition from linear per- 
turbations to the non-linear gravitational collapse of structures. The present 
catalogs do not yet see the transition since they do not look up sufficiently 
back in time. It can be noticed that some recent surveys begin to see a dif- 
ferent power-law behavior at large scales (A « 200 — AQQh~^ Mpc, e.g. Lin et 
al 1996, Scaramella et al. 1998). The interpretation depends however on the 
K-correction adopted, and the curvature of the Universe (Joyce et al. 1999). 

Summarizing, it appears that the high-mass cut-off could be at the present 
epoch at about 300 Mpc and 10^^ M©, the mass of the largest superclusters. 
The fractal extends over about 4 orders of magnitude is scale and 7 in mass. 
It is smaller than the ISM fractal, but since larger and larger structures could 
decouple from inflation and develop, this range could still increase. 



5.2 Correlation Function and Conditional Density 

The debate on the spatial extent of the fractal has been complexified by the use 
of the correlation function to quantify the scaling laws in galaxy distributions. 
The correlation function is defined as 

^, ^ _ < n{ri).n{ri + r) > ^ 
< n >2 

where n{r) is the number density of galaxies, and < ... > is the volume average 
(over d^ri). One can always define a correlation length ro by ^(rg) = 1. 

This definition involves the average density < n >, which depends on 
the scale for a fractal distribution. This is unfortunate, since the derived 
correlation parameters (slope and correlation lengths) then depend on the 
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galaxy sample used (see Davis 1997, Davis et al 1988, Pietroncro ct al 1997). 
The fact that there exists a correlation length does not mean that there is 
no fractal, because of its definition different from that used in physics; there 
^0 characterizes the exponential decay of correlations (^ e^^/'>") (for power 
decaying correlations, the correlation length is infinite). Davis & Peebles 
(1983) or Hamilton (1993) argue that the fractal of galaxies cannot have a 
large spatial extent, since the galaxy-galaxy correlation length ro is rather 
small. The most frequently reported value is rg ~ 5h~^ Mpc (where h = 
i/o/lOOkms-iMpc-i). 

The same problem occurs for the two-point correlation function of galaxy 
clusters; the corresponding ^(r) has the same power law as galaxies, their 
length To has been reported to be about rg « 25h~^ Mpc, and their cor- 
relation amplitude is therefore about 15 times higher than that of galaxies 
(Postman, Geller & Huchra 1986, Postman, Huchra & Geller 1992). The lat- 
ter is difficult to understand, unless there is a considerable difference between 
galaxies belonging to clusters and field galaxies (or morphological segrega- 
tion). The other obvious explanation is that the normalizing average density 
of the universe was then chosen lower. 

Assuming that the average density is a constant, while homogeneity is 
not yet reached, could perturb significantly the correlation function, and its 
slope, as shown by Coleman, Pietronero & Sanders (1988) and Coleman & 
Pietronero (1992). The function ^(r) has a power-law behaviour of slope —7 
for r < ro, then it turns down to zero rather quickly at the statitistical limit 
of the sample. This rapid fall leads to an over-estimate of the small-scale 7. 
Pietronero (1987) introduces the conditional density 

< n{r,).n{n+r) > 

r(r) = 

<n> 

which is the average density around an occupied point. For a fractal medium, 
where the mass depends on the size as M{r) oc r^, D being the fractal 
(Haussdorf) dimension, the conditional density behaves as r(r) ex r^^^. 
It is possible to retrieve the correlation function as 

< n > 
In the general use of ^(r), < n > is taken for a constant, and we can see that 

i:> = 3-7 . 

If for very small scales, both ^(r) and r(r) have the same power-law behaviour, 
with the same slope —7, then the slope appears to steepen for ^(r) when 
approaching the length tq. This explains why with a correct statistical analysis 
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(Di Nella et al 1996, Sylos Labini & Amendola 1996, Sylos Labini et al 1996), 
the actual 7 w 1 — 1.5 is smaller than that obtained using ^(r) (cf Fig. |^). 
This also explains why the amplitude of ^(r) and tq increases with the sample 
size, and for clusters as well. 

Note that the fractal distribution in the galaxy catalogs has been deter- 
mined from the light distribution, and this could be somewhat different from 
the mass distribution. There is also a morphological segregation of galaxies 
in clusters, and the different types of galaxies (ellipticals, spirals, or dwarfs) 
do not trace the same distribution. The mass distribution could be more 
complex. All this has led to the introduction of multifractality to represent 
the Universe (e.g. Sylos-Labini & Pietronero 1996). In a multifractal system, 
local scaling properties slightly evolve, and can be defined by a continuous 
distribution of exponents. This is a mere generalisation of a simple fractal, 
that links the space and mass distributions. Mutlifractality may also better 
account for the transition to homogeneity, with a fractal dimension varying 
with scale (Balian & Schaeffer 1989, Castagnoli & Provenzalc 1991, Martinez 
et al 1993, DubruUe & Lachieze-Rey 1994). 



5.3 Definition oj the Self- Gravity Domain 

The general view in cosmologiacl models is that the galaxy structures in the 
Universe have developped by gravitational collapse from primordial fluctua- 
tions. Once unstable, density fluctuations do not grow as fast as we are used 
to for Jeans instability (exponential), since they are slowed down by expan- 
sion. The rate of growth is instead a power-law, in the linear regime. If 5 is 
the density contrast: 

5{x) = (p(f)- < p >)/ < p > 

where < p > is the mean density of the Universe, assumed homogeneous at 
very large scale. If r is the physical coordinate, the comoving coordinate x is 
defined by r = a{t)x, where a(t) is the scale factor, accounting for the Hubble 
expansion (normalised to a(to) = 1 at the present time). Since the Hubble 
constant verifies H(t) = d/a, the peculiar velocity is defined by 

V ^ r — Hf = ax 

In comoving coordinates, the Poisson equation becomes: 

V^$ = 47rGa2(p- < p >) 

It can be shown easily that in a fiat universe, the density contrast in the 
linear regime grows as the scale factor a{t) = (t/to)^^'^- In the non-linear 
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Figure 5. Bottom: the average conditional density r(r) for several samples (Perseus-Pisces 
and CfAl in open rectangles, and LEDA in filled triangles, adapted from Sylos-Labini & 
Pictronero 1996). Top: ^(r) corresponding to the same surveys. The indicative line has a 
slope 7 = 1 (corresponding to a fractal dimension £> ~ 2). This shows that it is difficult to 
determine the slope on the ^(r) function. 



regime, on the contrary, the density p is much larger with respect to < p >, 
and the normal self-gravity and Jeans instability is retrieved. The receding 
velocities due to inflation do not exist anymore. The structures are bound 
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and decoupled from inflation. It is natural to define the self-gravity domain as 
limited by the largest decoupled structures, and if self-gravity is responsible 
for the fractal structure, it is expected that this scale will also delimit the 
transition to homogeneity. 



6 Statistical Mechanics of Self-Gravity 

Both for the ISM and for galaxy distributions in the Universe, self-similar 
structures are observed over large ranges in scales. Scaling laws are observed, 
which translate by an average density decreasing with scale as a power-law, of 
slope —7 between -1.5 and -1, corresponding to a fractal dimension _D = 3 — 7 
between 1.5 and 2. In the following, we will investigate the hypothesis that 
the main driver responsible for these fractal structures is self-gravity, even 
though the actual structures migh be more complex and perturbed. 

Since gravity is scale- independent, there are opportunities for a mecha- 
nism to propagate over scales in a self-similar fashion. For the ISM, in a quasi 
isothermal regime, a fractal structure could be build through recursive Jeans 
instability and fragmentation. This recursive fragmentation proceeds until 
the density is high enough to reach the adiabatic regime. Self-gravity could 
be the principal origin of the fractal, with generated turbulent motions in 
virial equilibrium at each scale. For galaxy formation, the smallest structures 
collapse first, and these influence the largest scale in a non-linear manner. It 
is obvious that in both cases, the system does not tend to a stationary point, 
but develops fluctuations at all scales, and these must be studied statistically. 

Recently, de Vega, Sanchez & Combes (1996a, b, 1998) have proposed a 
statistical field theory of self-gravity, not only to account for the existence 
of the structure, but also to be able to predict its fractal dimension and 
others critical exponents. This has been obtained by developping the grand 
partition function of the ensemble of self-gravitating particles. In transforming 
the partition function through a functional integral, it can be shown that 
the system is exactly equivalent to a scalar field theory. The theory does 
not diverge, since the system is considered only between two scale limits: 
the short-scale and large-scale cut-offs. Through a perturbative approach 
it can be demonstrated that the system has a critical behaviour, for any 
parameter (effective temperature and density). That is, we can consider the 
self-gravitating gaseous medium as correlated at any scale, as for the critical 
points phenomena in phase transitions (as was first suggested by Totsuji & 
Kihara 1969). 

Note that this approach is quite different from the thermodynamical ap- 
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proach developped by Saslaw & Hamilton (1984). Their theory is based on the 
thermodynamics of gravitating systems, which assumes quasi thermodynamic 
equihbrium. They justified this equihbrium at the smaU-scales of non-linear 
clustering, because the local relaxation and dynamical time-scales are much 
shorter than the expansion time-scale. The theory considers the essential pa- 
rameter b{t), the ratio of gravitational correlation energy to thermal kinetic 
energy, and deduces the value of this parameter from the observations. It ap- 
pears that b{t) varies also with scale. The predictions of the thermodynamical 
theory have been successfully compared with N-body simulations (Itoh et al. 
1993, Shcth & Saslaw 1996, Saslaw & Fang 1996). 



6.1 Hamiltonian of the Self- Gravitating Ensemble of N-bodies 

Let us consider a gas of particles submitted only to their self-gravity, in ther- 
mal equilibrium at temperature T {kT = /3^^). In the interstellar medium, 
quasi isothermality is justified, due to the very efficient cooling. For unper- 
turbed gas in the outer parts of galaxies, gas is in equilibrium with the cosmic 
background radiation at T w ZK (Pfenniger et al 1994, Pfenniger & Combes 
1994). For a system of collapsing structures in the universe, this can be a 
valid approximation, as soon as the gradient of temperature is small over a 
given scale. 

This isothermal character is essential for the description of the gravita- 
tional systems as critical systems, as will be shown later, so that the canonical 
ensemble appears the best adapted system. Moreover, the systems considered 
are not isolated gravitational systems, for which the microcanonical system 
should be used (e.g. Horwitz & Katz 1978a, b; Padmanabhan 1990). On the 
contrary, the mass or number of particles is not fixed, and according to the 
fluctuations, matter can enter any given scale, through condensation or evap- 
oration. The best statistical frame to consider is then the grand canonical 
ensemble, allowing for a variable number of particles N . The grand partition 
function Z and the Hamiltonian H^^ are 



oc /v „ „ W 



N\ J '" J ^^ h^ 

N=0 ■' ■' (=1 



<^^Vi d^li g-/3H^ 



Hn — 



N 2 1 

9m ^-^ \n. — 



l=\ l<l<j<N 



m 



m -<ij\ 
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where z is the fugacity — ea;p(— /3/ic) in terms of the gravito-chcmical potential 
^ic, and h is now the Planck constant. 

The latter integral can be transformed, using the continuous density 
p(^) " Si=i ^(f ^ Qj) and integrated to yield the potential, but intro- 
ducing a cutoff a for the minimum separation between particles, so that there 
is no problem of divergence. The cutoff a is here introduced naturally, it cor- 
responds to the size of the smallest fragments, or clumpuscules (of the order 
of ~ 10 AU). In fact, we consider that the particles of the system interact 
with the Newton law of gravity (l/r) only within the size range of the fractal, 
where self-gravity is predominant. At small scale, other forces enter into ac- 
count, and we can adopt a model of hard spheres to schematize them. Also at 
large scales, beyond the upper cutoff, different forces must be introduced. The 
phenomenological potential thus considered does not possess any singularity. 



6.2 Z as the Partition Function of a Single Scalar Field (j) 

Using the potential in l/r, and its inverse operator — ^V^ (but see also a 
similar derivation, with [1 — 0(a — r)]/r and its corresponding inverse operator, 
for the phenomenological potential with cutoff, in de Vega et al 1996b), the 
exponent of the potential energy can be represented as a functional integral 
(Stratonovich 1958, Hubbard 1959) 

IpGni-^ j :^^§^ p(x)p(y) ^ f f jy. -i^j d'xiWif +2,ny/^i j d'xi(x)p(x) 



With the change of variables: (f){x) = 2my/nGf3 ^(x) and 



/i3 -^^"-' ' -^J-r kT 

the partition function can be written as a functional integral Z = J J Dcp e^ 
where the action S is: 






1(V0)2 -^2 ^^(5) 



Note that the "equivalent" temperature Te// in the field theory is in fact 
inversely proportional to the physical temperature. It can be shown that the 
parameter /i is equal to the inverse of the Jeans length, itself of the order of 
the cutoff a. 

It is then possible to compute the statistical average of physical quanti- 
tites, such as the density p(r) 
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^eff J-eff 

where < . . . > means functional average over 0(.) with statistical weight 

eS[0(.)]. 

The equation for stationary points: V^cj) — —ji^ e"^*^^' has two main solu- 
tions. One is the constant stationary solution: (f>o — — oo; and the second is 
the singular isothermal sphere: (f>{r) = log — j-^- With a perturbative method, 
starting from the stationary solution cjjq — — oo, it has been shown that the 
theory scales, at large distances (de Vega et al 1996b). The same has been 
shown, starting from the isothermal sphere solution (Semelin et al. 1998). In 
the latter case, the development in series of the density correlation function 
is made with the relevant effective coupling constant: 

A = ^'^^ 



R 

and this coupling constant evolves through the rcnormalisation group equa- 
tions with the scale R, or scale ratio r — ln(R/a), where a is the low cut-off 
scale. A remarkable behaviour is found for A(t), since it vanishes periodi- 
cally, for values t„ = 27rn/\/7 {n integer). Periodically, the coupling constant 
diverges to infinity, which means that the perturbation theory is no longer 
valid, because of strong coupling. This occurs for scales: 

i?„ = i?oe2™/^ = i?o (10.749)" 

It appears then a hierachy of scale, varying by about an order of magnitude at 
each level. This numerical factor 10.749 depends essentially on the spherical 
geometry assumed for the computations, but is expected to be different for 
different geometries (like filaments, for example). 



7 Renormalization Group Methods 

The renormalization methods are very powerful to deal with self-similar sys- 
tems obeying scaling laws, like critical phenomena. In the latter case, exampli- 
fied by second order phase transitions, there exist critical divergences, where 
physical quantities become singular as power-laws of parameters called criti- 
cal exponents. These critical systems reveal a collective behaviour, organized 
from microscopic degrees of freedom, through giant fiuctuations and statistical 
correlations. Hierarchical structures are built up, coupling all scales together, 
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replacing an homogeneous system in a scale- invariant system. It can be shown 
that local forces are not important to describe the collective behaviour, which 
is only due to the statistical coupling of local interactions. Therefore, critical 
exponents depend only on the statistical distribution of microscopic configu- 
rations, i.e. on the dimensionalities or symmetries of the system. There exist 
wide universality classes, that allow to draw quantitative predictions on the 
system from only a qualitative knowledge of its properties (e.g. Parisi, 1988; 
Zinn-Justin 1989; Binney et al 1992). 



7.1 Critical Phenomena 

Critical phenomena occur at second order phase transitions, i.e. continuous 
transitions without latent heat. The paradigm of these systems is the transi- 
tion at the Curie point (T=Tc= 1043K) from paramagnetic iron, where the 
magnetic moment is proportional to the applied field m=/x B, to ferromag- 
netic state, where there exists a permanent magnetic moment iuq even in zero 
field. Although the permanent magnet tends to zero continuously at Tc, there 
are divergences: for instance the heat capacity C behaves as Coc \T — Tc|~", 
with a > 0. 

Also for the critical point of water, at which the transition from the liquid 
to gas becomes continuous, the compressibility kt cx: \T—Tc\~'^. At the critical 
point, it is easy to understand that the compressibility which tends to infinity 
generates large density fluctuations, and therefore light is strongly diffused by 
the varying optical index: this is the critical opalescence. The extraordinary 
fact is that microscopic forces can give rise to large-scale fluctuations, as if 
the medium was organized at all scales (cf Fig. 0). 

7.2 Universality of Critical Exponents 

Experiments have shown that the critical exponents for a wide variety of 
systems are the same, and more precisely they belong to universality classes, 
depending only on the dimensionality d of space and D of the order parameter 
(for instance if the fleld is scalar or a vector with dimension D). This univer- 
sality means that the details of the local forces are unimportant; therefore the 
local interactions can be simply modelled, through a schematic hamiltonian 
supposed to hold the relevant symmetries of the system. There exists only 2 
independent critical exponents. 

Renormalization techniques were applied to statistical physics in the 1970s 
(Wilson & Kogut 1974; Wilson 1975, 1983). In a renormalization group trans- 
formation, the scales are divided by a certain factor k, and particles are re- 
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Figure 6. Illustration of the fractal structure occuring when a Ising system passes through 
the critical temperature (a second order phase transition). The temperature is near zero at 
left, about critical in the center, and much larger than critical at right (from a simulation 
of the Ising model with a Monte-Carlo method) 



placed by a block of particles. Since the system is scale-independent, it should 
be possible to find an haniiltonian for the blocks which is of the same struc- 
ture as the original one. The new system is less critical than the previous one, 
since the correlation length S, has also been divided by k. It is a way to reduce 
the number of degrees of freedom of the system. With these techniques, pre- 
cise values of the critical exponents have been computed, for a whole range of 
models, and can be used for any problem in the same universality class. 



7.3 Statistical Self- Gravity 

As was shown in section 6, the grand-canonical self-gravitating system is crit- 
ical for a large range of the parameters, and it is difhcult to isolate a critical 
point, to identify diverging behaviours. However, it is well known (Wilson 
1975, Domb & Green 1976), that physical quantities diverge only for infinite 
volume systems, at the critical point. Since the self-gravitating systems are 
also finite and bounded, they only approach asymptotically the divergences. 
If A measures the distance to the critical point, (in spin systems for in- 
stance, A is proportional to \T — Tc\), the correlation length ^ diverges like 

e(A) ^ A--" 

and the specific heat (per unit volume) as C ^ A^" . But in fact, for a finite 
volume system, all physical quantities are finite at the critical point. When 
the typical size R of the system is large, the physical magnitudes take large 
values at the critical point, and the infinite volume theory is used to treat 
finite size systems at criticality. In particular, for our system, the correlation 
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length provides the relevant physical length ^ ^ R, and we can write 

The self-gravitating systems considered here have the symmetries d = 3 
and D = I (scalar field) , which should indicate the universality class to which 
it corresponds. It remains to identify the corresponding operators. Already 
in the previous sections, it was suggested that the field corresponds to the 
potential, and the mass density 

m p{x) = m e'^'^^'^ 

can be identified with the energy density in the renormalization group (also 
called the 'thermal perturbation operator'). 

We note that the state of zero density (or zero fugacity) , corresponds to 
a singular point, around which we develop the physical functions (and we 
choose A accordingly). At this point /i^/Te// = 0, the partition function Z is 
singular 

^2 /2nmkTY^^ 



i.e., the critical point A = corresponds to zero fugacity z. Writing Z as a 
function of the action S* at the critical point 



Z(A)= / f V^e-'''+^ I ""'"'* 



and computing statistical averages (de Vega et al. 1996a, b), the mass fiuctu- 
ations and corresponding dispersion can be found as: 

{AM{Rjf = <M^>-< M >2- / d^x d^y C{x, y) - R^'" 



AM{R) - R^/" 

This is the definition relation of the fractal, with dimension dn, and the 
scaling exponent v can be identified with the inverse Haussdorf dimension 
of the system, dn = -■ The velocity dispersion follows: Aw ^ i?'', with 

The scaling exponents i>, a have been computed through the renormal- 
ization group approach. The case of a single component (scalar) field has been 
extensively studied in the literature (Hasenfratz & Hasenfratz 1986, Morris 
1994a, b). For the Ising model d — 3, the exponent v = 0.631, from which 
we deduce dn = 1.585. Alternatively, in the case of weak perturbations, the 
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mean field theory can be applied, and dn — 2. Tliese values are compatible 
with the observed ones for astrophysical fractals. 



8 Conclusion 

We have emphasized the existence of two astrophysical fractals, the interstel- 
lar medium, with structures ranging from 10 AU to 100 pc, and the large-scale 
structures of galaxies, from 10 kpc to 200 Mpc at least. The first one is in 
statistical equilibrium, while the second one is still growing to larger scales. 
In both cases, we can describe these media as developping large-scale fluctu- 
ations with large correlations as is familiar in critical phenomena. We have 
investigated the hypothesis that in both cases, self-gravity is the main force 
governing these fractal structures. 

Numerical simulations can help to understand the formation of these 
structures, and the main mechanisms at play. Unfortunately, practical con- 
straints confine the possibilities to a limited range of scales, and results are 
often ambiguous. Simulations of MHD turbulence with self-gravity are not yet 
to the point to reproduce the scaling relations observed in the ISM. Purely self- 
gravitating simulations without large-scale injection of energy produce only 
transient fractal structures, depending in a large part on initial conditions. 
When large-scale energy is taken into account (by the galaxy shear namely), 
a quasi-stationary fractal structure, over ~ two orders of magnitudes, can be 
obtained. This is encouraging, waiting for more performant 3D simulations. 

A statistical thermodynamic approach of self-gravitating systems has been 
developped, and it is shown that the phenomenological potential, which is in 
l/r between two cutoffs (at small and large-scale), can be described by a 
scalar field theory. Using renormalization group methods, the system is found 
to be of the same universality class as the Ising d ~ 3 model. The critical 
exponents can then be derived, and the fractal dimension D — 1.6 deduced. 

The gravitational gas appears to be critical for a large range of temper- 
atures and couplings, while for spin models there is only a critical value of 
the temperature. This feature must be connected with the scale invariant 
character of the Newtonian force and its infinite range. 
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